M. W. Kim et al. 



Roles of orbital as a nexus between optical and magnetic properties in cubic i?Mn0 3 

(R = La, Pr, Nd, Gd, Tb) 

M. W. Kim 1 , S. J. Moon 1 , J. H. Jung 2 , Jaejun Yu 3 , Sachin Parashar 1 , P. Murugavel 1 , and T. W. Noh{] 

1 ReCOE & School of Physics, Seoul National University, Seoul 151-747, Korea 

2 Department of Physics, Inha University, Incheon 402-751, Korea 
3 CSCMR & School of Physics, Seoul National University, Seoul 151-747, Korea 

We investigated the afr-plane absorption spectra of i?Mn03 (R — La, Pr, Nd, Gd, and Tb) 
thin films. As the R-ion size decreases, we observed a drastic suppression of the 2 eV peak, i.e. 
the inter- site optical transition between spin- and orbital- aligned states across the Mott gap. We 
found that both lattice distortion and the corresponding orbital mixing of the ordered orbital state 
should play an important role in the 2 eV peak suppression. We also found that the 2 eV spectral 
weight is proportional to the A- type antiferromagnetic ordering temperature, which suggests that 
the magnetic interaction might be sensitively coupled to the orbital mixing. 

PACS numbers: 75.70.-i,77.90.+k,78.20.-e 



LaMnOs has been known as a mother compound of 
the colossal magnetoresistance(CMR) manganites, where 
charge, spin, lattice, and orbital degrees of freedom inter- 
play with each other to determine their intriguing physi- 
cal properties 0,0. LaMnOs has an orthorhombic struc- 
ture with four 3d-electrons: three t^g and one e g elec- 
trons. Since three t^g electrons form an orbitally closed 
shell, many physical properties are believed to be deter- 
mined by its e g electron. In its ground state, LaMnOs is 
a Mott-insulator 0, 0, with the A-type spin and the 
C-type orbital orderings, which are schematically drawn 
in Figs. 1(a) and 1(b), respectively. The antiferromag- 
netic (AFM) ordering temperature T/v is about 140 K, 
and the orbital ordering temperature is around 800 K. 
The occurrence of the spin- and orbital-ordered state has 
been understood in terms of the cooperative Jahn- Teller 
(JT) transition 0. 

Rare-earth substitutions of the La ion provide an in- 
triguing phase diagram for i?Mn03 (R = rare-earth ion) 
0, as illustrated in Fig. 1(c). As the i?-ion size tr 
decreases, the crystal structure of RM11O3 changes from 
orthorhombic (R = La-Dy) to hexagonal (R = Ho-Lu). 
TbMnOs and DyMnOs are located near the structural 
phase boundary, and they have attracted lots of atten- 
tion recently due to their complicated low temperature 
magnetic states and multiferroic properties 8]. On the 
other hand, the magnetic properties of the orthorhom- 
bic perovskite i?Mn03 (R = La-Tb) has a rather simple 
independence: Tn decreases with decreasing tr. Struc- 
tural deformations, such as buckling and distortion of 
the Mn06 octahedra, also increase. According to the 
Goodenough-Kanamori rule 0, , the orbital overlap of 
electrons should be crucial in determining the magnetic 
interaction. The rule takes into account the overlap in 
terms of the Mn-O-Mn bond angle <fi. However, the rapid 
decrease of T/v with the R-ion substitution is rather unex- 
pected, since the change of <p is less than 10°. In addition 
to the lattice distortion due to the <p variation, a neutron 
scattering measurement showed that the e g electron state 



should have a mixed character of \3z 2 — r 2 ) and \x 2 — y 2 ) 
orbitals and that the degree of orbital mixing varies sys- 
tematically with tr This orbital mixing could also 
affect the orbital overlap of electrons. Therefore, RMnOs 
(R = La-Tb) is an ideal system to investigate roles of 
the lattice distortion and the orbital mixing in numerous 
physical properties. 

Optical spectroscopy has been known to be a pow- 
erful tool to investigate the orbital degrees of freedom 
001 EH [13 In this Letter, we report the afr-plane opti- 
cal responses of epitaxial i?Mn03 (R = La, Pr, Nd, Gd, 
and Tb) films. We find that the spectral weight of the op- 
tical transition across the Mott gap, located around 2 eV, 
decreases rapidly as tr decreases. This dramatic reduc- 
tion of the spectral weight cannot be explained in terms 
of the conventional model based on the structural varia- 
tions. We demonstrate that the spectral weight change 
could be explained by taking account of the structural 
distortion, i.e., the change of 0, and the orbital mixing. 
We also find that the measured spectral weight change is 
proportional to the variation of T/v. 

High quality i?Mn0 3 (i?=La, Pr, Nd, Gd, and 
Tb) thin films were grown on double-side-polished 
(LaA103)o.3(SrAlo.5Tao.503)o.7 substrates by using the 
pulsed laser deposition. From x-ray diffraction measure- 
ments, it was found that all the films grew epitaxially 
with their c-axis perpendicular to the film surfaces. De- 
tails of the film growth and their characterization were 
reported elsewhere [l3| . Transmission spectra of the films 
were measured from 0.4 to 4.0 eV by using a grating 
spectrophotometer. The absorption coefficients were de- 
termined by taking the logarithm of the transmittance, 
subtracting that of the substrate, and dividing by the 
film thickness. Since the normal-incident optical geom- 
etry was used, the absorption spectra should come from 
the a6-plane responses of the films. Figure 2(a) shows 
the absorption spectra ql(oj) of RMnO^ at room tem- 
perature, where all the samples should be in the C-type 
orbital ordered state 14]. The spectra of LaMnOs are 
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FIG. 1: (color online), (a) The A- type spin and (b) the 
C-type orbital ordering pattern of electrons at the Mn-sites 
in LaMnOs. (c) A schematic magnetic phase diagram of 
i?Mn03, redrawn from Ref. 7]. The spectral weight of 2 eV 
peak, Sw (the solid square), shows similar R-ion size depen- 
dence with the A-type AFM ordering temperature (T) (the 
open circles). The i^-type AFM ordering T (the open dia- 
monds) and the incommensurate spin ordering T (the open 
triangles) are also shown. The hatched area represents the 
region where the commensurate spin order and ferroelectric 
property emerges. 



composed of a peak near 2 eV and much stronger ab- 
sorption peaks above 3 eV. The higher energy absorp- 
tion features come from the charge transfer transition 
from O 2p to Mn 3d 05 El- After long debates, numer- 
ous recent experiments clearly demonstrated that the 2 
eV peak should be interpreted as an inter-site transition 
across the Mott gap in the orbitally degenerate Hub- 
bard model (ODHM) 0, H H As shown in Figs. 1(a) 
and 1(b), the correlation-induced transition within the 
afr-plane should occur between the e g electron states at 
the neighboring sites in the ferromagnetic-s pin (FM) and 
antiferro-orbital (AFO) configuration kJ&Oil The ab- 
sorption spectra of other i?Mn03 have very similar spec- 
tral features. As tr decreases, a(uo) for the charge trans- 
fer transition above 3 eV is nearly independent of the R- 
ion, however, ql(oj) for the correlation-induced 2 eV peak 



becomes strongly suppressed. To obtain more quantita- 
tive information, we estimated the spectral weight Sw 
by subtracting the charge transfer transition background 
and integrating a(co) from 0.2 to 2.7 eV. The experimen- 
tally determined Sw, marked as the solid squares in Fig. 
2(b), becomes drastically suppressed with decreasing vr. 

Such a dramatic decrease of Sw is rather unexpected. 
All the RM11O3 compounds, studied in this work, have 
the same orthorhombic crystal structure and the same 
spin/orbital ordering pattern, but only with a relatively 
small variation of <p. Let us look into the possible role 
of the structural variations of RMnOs in the large Sw 
change. According to the chemical grip estimate [l7|, 
the inter-site transition between the d states can vary 
approximately as cos A (j). As shown in Fig. 2(b), the con- 
tribution of the structural variations to the Sw change 
could be as large as 30%, but is still much smaller than 
the experimentally observed Sw changes. Therefore, 
the electronic structure change due to the structural 
variation alone cannot explain the large suppression of 
Sv^. To elucidate the origin of the 2 eV peak spectral 
change, we applied the Fermi- Golden rule and evaluate 
the corresponding matrix element. Within the electric 
dipole approximation Sw becomes proportional to 
\(ipf\ — ih\7\i/ji)\ 2 . This matrix element can be approxi- 
mated by the second order perturbation, similarly to the 
superexchange process, through the oxygen p orbit als: 

s w cx K^/lviVi)! 2 ~ l£»a><Pak>l7A, (i) 



by assuming that the energy gap A remains almost 
unchanged [13 • Here |^) and \i/jf) represent the wave- 
functions of the initial occupied and the final unoccupied 
Mn e g orbitals, respectively, of the transition considered. 
\Pa) {ol = x, y, z) represents the oxygen p orbitals which 
bridge the Mn e g orbitals. For the 2 eV peak, the ma- 
trix element can be estimated by the inter-site transition 
from the occupied e g state at one site to the unoccu- 
pied e g state at the neighboring site with the FM/AFO 
configuration, as shown in the inset of Fig. 2(a). 

First, let us consider the contribution of the Mn-O- 
Mn bond angle change in the matrix element. As shown 
in Fig. 3(a), the buckling of the Mn06 octahedra in the 
GdFeOs type lattice will cause a decrease in (j). With- 
out the buckling {i.e. (j) ~ 180°), the larger lobe of the 
\3x\ — r\ )-type orbital of a Mn 3+ site is aligned to face 
the smaller lobe of the neighboring orbital orthogonally, 
as shown in Fig. 1(b). As the buckling is turned on, the 
orbital lobe of an electron at one site will rotate with 
respect to that at the neighboring site, which will result 
in a reduction in the inter-site hopping amplitude and 
thereby a decrease in S^. To evaluate the changes in 
|^.) and 1^) quantitatively, the rotation of the orbitals 
in the afr-plane was formulated in terms of the rotational 
transformation of the local (#2, V2) coordinates by (f) with 
respect to the local y{) coordinates, as shown in Fig. 
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FIG. 2: (a) Absorption spectra of i?Mn0 3 (R = La, Pr, Nd, 
Gd, and Tb) thin films at room T. The inset schematically 
represents the inter-site transition corresponding to the 2 eV 
peak, (b) The experimental Sv^ (the solid squares) is com- 
pared with the calculation results: the simple estimation of 
bandwidth change (~ cos 4 (f)) (the open circles), our model cal- 
culation for various (j) by using the orbitals of LaMnOs, i-e., 
at a fixed (= 107.6°) (the open triangles), and our model 
calculation for various and (the open diamonds). 



3(a). Here, the local z\ and z^ axis directions are as- 
sumed to be the same. When the orbital wavefunctions 
of LaMnOs were used, it was found that the orbital rota- 
tion effect on Sw is proportional to cos 2 (tt— <j)). As shown 
in Fig. 1(c), the R-\on substitution in i^MnOs makes <p 
vary from 155.2° to 145.3° [Ullj. In Fig. 2(b), the cal- 
culated values of Sw are plotted with the open triangles. 
It is obvious that the variation in <ft alone cannot account 
for the large change in the experimental Sw- 

Now, let us include the orbital mixing contribution. In 
a cubic Mn06 octahedron, two e g orbitals remain dou- 
bly degenerate. Under the JT-type distortion along the 
z-direction, the e g orbitals become split into two orthog- 
onal orbitals, i.e. \3z 2 — r 2 ) and \x 2 — y 2 ). However, 
a neutron scattering experiment showed that the actual 
occupied e g orbital of RM11O3 should be a mixed state 
of these two orbitals depending on the local distortion of 
the Mn06 octahedron, and further that the degree of the 



orbital mixing will vary depending on tr To include 
the orbital mixing effects in Eq. (1), we constructed re- 
alistic Mn e g orbitals using the orbital mixing angle 0: 
the occupied orbital at site 1 and the unoccupied orbital 
at the neighboring site 2 are written as 
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where the subscripts in the wavefunctions represent the 
different local coordinates. The unoccupied orbital, cor- 
responding to the final state of the transition, is orthog- 
onal to the occupied orbital at site 2. To visualize the 
orbital mixing effects, the occupied orbitals for three dif- 
ferent values are plotted in Fig. 3(b). [Note that the 
orbital configuration shown in Fig. 1(b) corresponds to 
<j> = 180°, = 108°.] From Eqs. (1) and (2), we obtained, 
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Using the reported (0>4>) values from the neutron scat- 
tering experiment [l(J, l2i| ■> we can estimate the values of 
Syr (#,(/>) and plot them with the open diamonds in Fig. 
2(b). The estimated (#,(/>) values agree quite well with 
the measured Sw change, indicating the importance of 
the orbital mixing. The 0- and (^-dependence of is 
displayed in Fig. 3(c). Note that = 107.6° for LaMn0 3 
and = 114.3° for TbMn0 3 . Although the variation 
of the value is about 6.7°, smaller than that of the <j) 
value (i.e. about 9.9°), the variation of Sw due to the 
change in the orbital mixing is larger than that due to 
the (f) change. A possible reason is the strong anisotropy 
in the shape of the orbitals. When 9 = 90°, the oc- 
cupied and unoccupied orbitals given in Eq. (3) have 
the mean state of the two orthogonal orbitals. As 6 in- 
creases, the \x 2 — y 2 ) orbital enhances the wavefunction 
overlap and the \3z 2 — r 2 ) orbital reduces it within the 
afr-plane, so the mixing of those two orbitals results in 
the minimum around = 120° in Eq. (3). As shown in 
Fig. 3(c), i^MnOs are located near the (0,0) space where 
Sw will change rapidly and depend strongly on 0. Thus, 
the orbital mixing becomes a crucial factor in numerous 
physical properties of RMxiOz, including the change in 
Sw- 

The Sw values for various manganites are marked with 
the solid squares in Fig. 1(c). It is remarkable to note 
that the ^-dependence of the Sw change is quite simi- 
lar to that for T/v, the A- type AFM ordering tem- 
perature. With decreasing R-ion radius (tr) from La 
to Tb, T/v systematically decreases, i.e. from ~ 140 
K for LaMnOs to nearly zero for TbMnOs. Similarly, 
Sw also becomes significantly reduced for TbMnOs as 
compared with the value for LaMnOs. These similar in- 
dependences of Sw and T/v are rather surprising, since 
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FIG. 3: (color online). The setup is schematically shown for 
the calculation of the orbital mixing angle (9) and the Mn-O- 
Mn bond angle (0) dependent Sv^. (a) The local coordinates 
(x,y) for the Mn06 octahedra and the Mn e g orbitals (The 
orbital lobes drawn with the dashed lines represent the un- 
occupied orbitals.) (b) The occupied orbitals of selected 
values, (c) The calculated Sw as a function of and <fi. 



the Sw change comes from the ab-plane response while 
the AFM ordering at T/v occurs along the c-axis. 

One possible explanation for this intriguing phe- 
nomenon could be an occurrence of additional FM com- 
ponent in the inter-plane interactions due to the buckling 
of the Mn06 octahedra. In the undistorted case (i.e., 
<j) = 180°), the A-type spin order should occur due to the 
FM e g -e g interaction within the afr-plane and the AFM 
t2g-t2g superexchange interaction along the c-axis. Here 
T/v is mainly determined by the latter, since it is much 
weaker than the afr-plane FM interaction [2]] . When the 
buckling of the Mn06 octahedra occurs, however, the 
overlap of the e^-orbitals between the Mn-planes brings 
out a new FM interaction along the c-axis. Although 
Tat of i?Mn03 should be determined by the competi- 
tion between the AFM and the FM interactions, the tr 
dependence of T/v could be realized mostly by the lat- 
ter: the AFM interaction should not be so sensitive to 
the structural change due to the nature of the orbitally 



closed t2 9 levels, but the FM interaction should critically 
depend on the buckling. Since the FM interaction does 
appear from the tilting of the Mn06 octahedra, it could 
be closely related to S^, which will be proportional to 
the square of the electron hopping matrix in the a6-plane. 
This scenario suggests that the AFM and FM interac- 
tions will compete with each other and achieve a bal- 
ance around TbMnOs. To explain the anomalous mag- 
netic ground states near the phase boundary, shown in 
Fig. 1(c), Kimura et al. recently used a two-dimensional 
anisotropic neighbor interaction model [2^. Our picture 
based on a new FM interaction along the c-axis might 
provide an alternative starting point to explain the in- 
triguing magnetic states near the multiferroic phases. 

The conventional Goodenough-Kanamori rule takes 
into account of the orbital overlap in terms of 0, El • 
Then, the sign of the effective magnetic interactions {i.e. 
AFM and FM ground states) is expected to change near 
(j) = 135°. In i?Mn03, the GdFeOs type distortion can 
induce the competition between AFM and FM along the 
c-axis, so one could envisage a disappearance of AFM 
near 135°. However, as displayed in Fig. 1 (c), the A- 
type AFM order in the orthorhombic disappears 
at ^ 145°. This deviation from the Goodenough- 
Kanamori rule should originate from the additional con- 
tribution of the orbital mixing to the orbital overlap. 

In summary, we reported that the spectral weight of 
the 2 eV peak changes drastically with rare-earth ion size 
in i?Mn0 3 (R=La, Pr, Nd, Gd, and Tb). The spectral 
weight change was successfully understood in terms of 
the optical matrix element in which Jahn- Teller distor- 
tion and the rotation of the orbital were taken into ac- 
count within the orbitally degenerate Hubbard picture. 
Similar behaviors between the 2 eV spectral weight and 
the A-type antiferromagnetic ordering temperature sug- 
gest that the superexchange interaction in i?Mn03 might 
be tuned by the orbital degree of freedom. 
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